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A new method for estimating frequency-dependent core 
shifts in AGN jets. 
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ABSTRACT 

We discuss the opacity in the core regions of Active Galactic Nuclei (AGN) observed 
with Very Long Baseline radio Interferometry (VLBI) and describe a new method 
for deriving the frequency-dependent shifts of the VLBI core from the frequency- 
dependent time lags of flares observed with single-dish observations. Application of 
the method to the core shifts of the quasar 3C 345 shows a very good agreement 
between the core shifts directly measured from VLBI observations and derived from 
flares in the the total flux-density using the proposed method. This provides direct 
evidence that the observed jet component speeds in this AGN represent the actual 
physical speed of the feature, rather than a pattern speed. The frequency-dependent 
time lags of flares can be used to derive physical parameters of the jets, such as distance 
from the VLBI core to the base of the jet and the magnetic fields in the core region. 
Our estimates for 3C 345 indicate core magnetic fields ~ 0.1 G and magnetic field at 
lpc~0.4G. 
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1 INTRODUCTION 

Most Active Galactic Nuclei (AGN) display a "core+jet" 
structure in Very Long Baseline Interferometry (VLBI) im- 
ages. In the standard interpretation, the core is taken to 
be th e optically thick base of the jet (Bl andford fc Konigll 
1 19791 ). Due to synchrotron self- absorption, the absolute po- 
sition of the observed VLBI core (r = 1 surface) shifts 
systematically with frequency, moving in creasingly ou tward 
along the VLBI jet with lower frequency (Konigl 1981). This 
frequency-dependent core shift has a direct effect on astro- 
metric measurements performed in the radio and optical. 
In the near future, the GAIA astrometry mission and the 
Space Interferometry Mission will begin. For both missions, 
matching the o ptical astromet ric catalogues to the radio 
catalogues fe.g. lFev et ai1l200ll ) presents a very important 
problem. The core shifts can introduce offsets between the 
optical a nd radio positions of AGN of up to several milliarc- 
seconds (|Kovalev et al.ll200Si ). which will strongly affect the 
accuracy of matching the radio and optical catalogues. 

Core shifts are also needed for the correct reconstruction 
of VLBI spectral-index and rotation-measure maps. More- 
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over, knowledge of the frequency-dependent core shifts can 
be used to derive physical parameters of the jet, such as the 
core magnetic field and the distance from the VL BI core to 
the base of the jet |Lobanovlll99Sl ; lHirotani|[2005l ). 

Therefore, precise measurements of the core shifts 
are necessary. One way to obtain core shifts is through 
phase-referencing VLBI observations, however this is a 
complex and resource-intensive technique, and phase- 
referencing core shifts have b een determined for only a 
few AGN, s uch as 1038+5 2 8 (iMarcaide fc Shapiro! \l9m), 
4C 39.25 l|Guirado et all Il995l). 3C 39 5 dLara et all 
1 19941) 3C 390.1 jRos fc Lobanovl l200ll ). and 
M 81 (|Bietenholz, Bartel. fc RupenlbOoi T Another indirect 
method to measure core shifts is to align opti cally thin parts 



of an AGN jet at differen t frequ en cies (e.g. iKovalev et al 
120081 ; ICroke fc Gabuzdal 120081 ; lO'Sullivan fc Gabuzda 
120091) . However, this method requires simultaneous multi- 
frequency VLBI observations, which are likewise fairly 
resource intensive, and does not always yield unambiguous 
results. The limitations of these techniques are exacerbated 
by the fact that the core shift may well depend on the 
activity state of an AGN, and therefore be time dependent, 
whereas at present only isolated core-shift measurements for 
individual AGN are available. In addition, the magnitude 
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of the core shifts that can be detected is limited by the 
resolution of the VLBI observations used. 

In this paper, we discuss the opacity in the core regions 
of AGN and describe a new method for deriving core shifts 
from the frequency-dependent time lags of flares observed 
with single-dish observations. We use these time lags to de- 
rive physical parameters of the jets, such as the distance 
from the VLBI core to the base of the jet and the magnetic 
fields in the core region. 



2 FREQUENCY-DEPENDENT TIME LAGS 

The core is a compact feature in the VLBI map of an AGN 
with high brightness and a relatively flat spectrum, which is 
usually interpreted as that part of the jet where the optical 
depth is r = 1. Since the r = 1 surface h as different loca- 
tions at different frequencies (Konigl 198l]), the absolute po- 
sition of the core should depend on the observing frequency. 
According to the standard shock-in-jet model, a change of 
electron density and pressure at the injection point near the 
base of the jet will cause a shock wave to propagate along 
the jet. This will cause brightening of the core region, which 
will be seen as a flare in the total flux-density light curves, 
and will be followed by the appear ance of a jet component 
(or components) in VLBI maps (e.g. iMarscher fc Gearil 19851 ; 
I Gomez et al.l fl997l ). Consider a conical jet geometry (see 
Figure [1} observed at a viewing angle (p. A shock wave ap- 
pears at the distance R on and then propagates along the jet. 
It will cross the r = 1 surface for a particular frequency Vi 
at some time T% and emerge out of the core region. Since 
the position of the r = 1 surface is shifted further from the 
jet base at lower frequencies, the times T< will be delayed 
at lower frequencies. Figure [1] shows how the T; correspond 
to the maxima of the total flux-density flares obtained with 
single-dish observations. Crossing the r = 1 surface corre- 
sponds to a maximum of the total flux-density outburst at 
the corresponding frequency, and the time of the flare max- 
imum thus also depends on the frequency; thus, the time 
lags contain information about core opacity. According to 
the geometry of Figure [T] and the principles of superlumi- 
nal m otion (e.g. lReeslll967l ; pTurler. Courvoisier. fc Paltanil 
2000), the time t that has passed after the shock wave's 
appearance at the distance R on in the observer's frame is 



t 



(1 + z)simp{R(v) — R on (u)) 

PappC- 



(i) 



where 7? is the distance along the jet axis in the rest frame 
of the quasar, /3 app is the apparent velocity in the plane of 
the sky in units of c, <p is the viewing angle, and z is the 
redshift of the source. In this equation, we assume that the 
apparent speed is the actual speed of the shock. Observing 
at multiple frequencies v a and v b , we can estimate the time 
lags between the maxima of the total flux-density outburst: 



At 



t(v a ) - t{Vb) 
(1 + z)sinip 



(2) 



[R(v a ) ~ R{fb) + Ron(v a ) - R on (vd) 



Assuming that the shock wave appears at the same distance 
-Ron at all frequencies, we obtain At(v) obs — AR pro j / Pa PP c, 
where AR pro j is the distance between R at the two frequen- 
cies projected onto the plane of the sky. Taking expression 



33) f or the frequency-dependent core shift from iHirotanil 
2005), we can write the dependence of the time lags on 

the frequency, spectral index a, magnetic field and electron 

density: 



At(u) oba 



sirup t k 

Pa 



l/*v 



l/k r 



l/k r 



l/k r l/kr 



where 



(3 - 2a)m + 2n - 2 



2a 



3 - 2a 
5 - 2a 



(4) 



(•>) 



(6) 



xb is a dimensionless variable for B and / is a function of 
N\, the spectral index, and the viewing angle. 

Defining the core-position (or time-lag) offset as 



Q, rv = 4.85 x 10" 



i At/3 app DL 
(1 + z) 2 



l/k r l/kr 



l/kr 



l/kr 



(l) 



we can obtain formulas for the magnetic field and the dis- 
tance from the VLBI core to the central engine. In this for- 
mula, Q rl , is measured in units of pc ■ GHz, At in yrs, j3 apP 
in milliarcsecond (mas) per year, and Dl is the luminosity 
distance of the source in parsec. 

Formula [3] assumes that the electron number density 
and magnetic field scale with distance along the jet R as 



N* = N ir - n ,B = B ir ~ 



(8) 



where Ni and Bi refer to the v alues of N* and B at r — 1 pc. 

Following iLobanovl (|l998l ) and IHirotanil l|2005l ). we can 
write the distance between the core and the base of the jet 
in terms of the frequency-dependent time lags as 



/ \ Qrv —l/k r 
r C ore\V) = —. V , 

sirup 



(9) 



where k r is estimated using the frequency-dependent time 
lags. We can calculate the magnetic field at 1 pc distance 
for the special case when there is equipartition between the 
energies of the particles an d magnetic field (k r = 1) and the 
spectral index is a = -0.5 l|Q'Sullivan fc Gabuzdall2009h : 



Bi = 0.025 



n^ji + z) 2 

S 2 8sin 2 ip 



1/4 



(10) 



where 5 is the Doppler factor, 9 the jet half-opening angle, 
ip the viewing angle, and Bi is in Gauss. The equipartition 
magnetic-field strength in the core can then be found from 
the relation 

B core {v) = BirZove (11) 
for B core at a particular frequency v. 



3 FITTING THE FLUX-DENSITY LIGHT 
CURVES 

In order to check the proposed method for measuring the 
frequency-dependent core shifts from the time lags and the 
validity of our assumptions, we calculated time lags from 
the total flux-density light curves for the AGN 3C 345, 
whose core shift has been measured at the same frequencies 
used to construct the light curves. The frequency-dependent 
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Figure 1. Sketch of the jet geometry discussed in section [2] as observed in the rest frame of the quasar. The shaded rectangle shows a 
shock wave that starts at the distance R on and propagates along the jet axis with the speed j3. The vertical lines mark the core positions 
at different frequencies and the times when the shock wave reaches the core at a particular frequency. The times Tj correspond to the 
maxima of the flares at frequencies i/;. 
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Figure 2. Total flux-density light curve of 3C 345 at 4.8 GHz, 8 GHz, 14.5 GHz, 22 GHz, and 37 GHz. The data are from the University of 
Michigan and Mctsahovi Radio Observatory monitoring databases. The arrows mark the times of the ejections of the new jet components 
C2, C4 and C7. 
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Figure 3. Total flux-density light curve of 3C 345 at 14.5 GHz 
(the data are from the University of Michigan monitoring 
database). The long dashed lines show decomposed Gaussian 
components and the solid line the sum of the fitted Gaussians. 



time delays were calculated by fitting Gaussian functions 
to the total flux-density light curves from the University 



to tne total flux-densit y ugnt cu rves tro m tne Um\ 
of Michigan (|Aller et al.ll 19991 ; lAller Aller. fc Hughes 



and Metsahovi R adio Observatory ( Terasranta et al 



2003) 



1992 



Il998l , l2004l , l2005l ) monitoring databases at 4.8 GHz, 8 GHz, 
14.5 GHz, 22 GHz, and 37 GHz. 

The total flux-density light curves were decomposed 
int o Gaussian compo n ents f ollowin g the procedure d is cusse d 
in iPvatunina et al.l (|2006l ) and IPvatunina et all <|2007l ). 
Long-term trends in the total flux-density variations were 
calculated as polynomial approximations for the deepest 
minima in the light curves and subtracted before the fit- 
ting. The Gaussian decomposition was performed such that 
it first removed the trend, then found the highest peak in 
the light curve and fitted a Gaussian to the peak, based on 
a x 2 minimization. The fitted component was then removed 
from the light curve and the procedure was repeated until 
all significant peaks were fitted with Gaussians. During the 
fitting, we applied the general criterion that the smallest 
number of individual flares (Gaussian components) provid- 
ing a complete description of the light curve was used. The 
number of fitted Gaussians depended on the time interval 
covered by the light curve and the characteristic time scale 
of the source variability. 



4 TEST OF THE METHOD: THE QUASAR 
3C 345 



The quasar 3C 345 (z = 0.5928, iMarziani et~aH (|l996l )) 
is one of the best studied AGN on VLBI scales. Sev- 
eral jet components with apparent velocities of 2 — 20c 
have bee n observed, moving along strongly curved trajecto- 
ries (e.g. lUnwin et al.lll983l: IZensus. Cohen, fc Unwinlll995l: 
Lobanov fc Zensusl Il999l; iRantakvro. Baath. fc Matveenkd 
19951 ; iRos. Zensus. fc Lobanovll2000l ) ? found evidence for 
periodic changes in the parsec-scale jet with a 9-year period. 

The total flux-density light curves at 4.8 GHz, 8 GHz, 
14.5 GHz, 22 GHz and 37 GHz are shown in Fig. H The 



best x 2 values were reached with a fit of 9 Gaussians to the 
light curve; the Gaussians fitted are shown together with the 
light curve at 14.5 GHz in Fig. [3] The long-dashed curves 
show the Gaussians and the solid curve the final sum of all 
9 fitted Gaussians. It is clear that the sum of the Gaussians 
fits well all the main features in the light curve. The param- 
eters of the 9 Gaussians are given in Table [1] The columns 
of this table give the (1) component designation, (2) observ- 
ing frequency, (3) epoch of maximum flux, (4) maximum 
amplitude, (5) full width at half maximum (FWHM) of the 
Gaussian component corresponding to the outburst O, (6) 
time delay between the epoch of maximum flux for the out- 
bursts at the given frequency and the corresponding epoch at 
the highest available frequency AT, (7) calculated spectral 
indexes and (8) calculated k r values (the k r values will be 
discussed below). We use a positive spectral index conven- 
tion S ~ v a . We measured the spectral index for each flare, 
fitting a linear regression into log(f) vs. log(S) plots, where 
S is an amplitude of a flare, obtained as the maximum of 
fitted Gaussian function. Spectral index was calculated only 
for those outbursts that have been detected at three or more 
frequencies. 

We were able to calculate time-delay core shifts for the 
epochs of three outbursts associated with the ejection of 
components for which superluminal speeds have been pub- 
lished: C2, C4 and C7, marked in Fig. Core-shift mea- 
surements based on a direct comparison of VLBI images at 
different frequencies are available for epochs near our mea- 
sured time delays. Table f2] summarizes and compares the 
core shifts directly measured from the VLBI images and our 
calculated core shifts derived from time lags obtained with 
single dish observations. The epoch of the VLBI observa- 
tions and name of the ejected jet components used for the 
ana lysis are listed in this table as we ll. 

iBiretta. Moore. &: Cohenl l| 19861 ) measured core shifts 
for 3C 345 of 0.05±0.13 mas between 10.7 GHz and 5 GHz in 
1982 and 0.111 ±0.007 mas between 22.2 GHz and 10.7 GHz 
in 1983.5, during two powerful flares that were associated 
with the ejection of the new jet components C2 and C4 

(Fig. ED. 

Compone nt C2 displayed an apparent speed of 0.48 ± 
0.02 mas/yr (|Biretta et al.lll986 | ); similar speeds were de- 
tected for this component by Zensus et alj |l995l ). fi app = 
0.4-0.53 mas/yr. The measured time lag between 5 GHz 
and 8 GHz for the 1982 outburst corresponding to the ejec- 
tion of C2 is 0.12 ± 0.04 yrs. Following our method of cal- 
culating the core shifts from the time lags as AR pro j [mas] 
= At(v) b s [yrs] -Papp [mas/yr], this corresponds to a core 
shift of 0.06 ± 0.02 mas. This "time-lag" core shift coincides 
well with the core shift of 0.05 ± 0.13 mea sured directly by 
aligning VLBI images (|Biretta et al.lll986l ). 

The jet component C4 displayed an appare n t spee d 
of 0.225 ± 0-015 mas /vr (ICaproni fc Abraham! l2004al ). 
I Zensus et alj (| 19951 ) and lBiretta et al.l (| 19861 ) detected simi- 
lar speeds of 0.32±0.15 mas/yr and 0.295±0.009 mas/yr for 
this jet component. The calculated time lag between 22 GHz 
and 8 GHz for the 1984 flare is 0.45 ± 0.08 yrs. If we use an 
average of the three speed estimates, this yields a time-delay 
core shift of 0.12±0.02, which is very close to the core shi ft 
of .111 ± 0.007 dire ctly measured bv lBiretta et al.1 (|l986l ). 

iLobanovl (|l998h measured core shifts of Ar = 0.05 ± 
0.03 mas (5.0-8.4 GHz), Ar = 0.21 ± 0.06 mas (8.4 - 
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Table 1. 3C 345: Parameters of outbursts 



Comp. 


Frcq. 
(GHz) 


Amplitude 

(Jy) 


Tvnax 

(yr) 




(yr) 


Time delay 

(yr) 


Sp. Index 


k r 


A 


37.0 


12.69 ± 0.10 


1982.19 ± 0.02 


3.90 ± 0.05 


0.00 ± 0.04 


0.30 ± 0.07 


1.75 ± 0.02 




22.0 


11.79 ± 0.06 


1982.29 ± 0.07 


6.05 ± 0.08 


0.11 ± 0.09 








14.5 


11.35 ± 0.04 


1982.27 ± 0.01 


3.87 ± 0.02 


0.09 ± 0.03 








8.0 


9.49 ± 0.01 


1982.74 ± 0.01 


4.81 ± 0.01 


0.55 ± 0.03 








4.8 


6.56 ± 0.03 


1982.86 ± 0.03 


3.92 ± 0.03 


0.67 ± 0.05 






B 


14.5 


5.71 ± 0.03 


1985.19 ± 0.01 


1.86 ± 0.01 


0.00 ± 0.03 


0.98 ± 0.79 






8.0 


1.52 ± 0.02 


1985.05 ± 0.01 


1.25 ± 0.01 


0.14 ±0.02 








4.8 


2.01 ± 0.02 


1985.57 ± 0.01 


1.18 ± 0.01 


0.38 ±0.02 






C 


37.0 


5.54 ± 0.04 


1986.58 ± 0.01 


2.53 ± 0.03 


0.00 ± 0.02 




1.81 ± 0.05 




22.0 


3.45 ± 0.04 


1987.05 ± 0.04 


3.38 ± 0.04 


0.47 ± 0.05 








14.5 


5.88 ± 0.03 


1987.07 ± 0.01 


1.84 ± 0.01 


0.49 ± 0.02 








8.0 


5.12 ± 0.02 


1987.09 ±0.01 


3.62 ± 0.02 


0.51 ± 0.02 








4.80 


5.39 ± 0.01 


1987.38 ± 0.01 


4.10 ± 0.01 


0.81 ± 0.02 






D 


37.0 


1.56 ± 0.04 


1988.70 ± 0.01 


0.90 ± 0.01 


0.00 ± 0.02 








14.5 


1.70 ± 0.02 


1988.80 ± 0.01 


1.17 ± 0.01 


0.10 ± 0.02 






E 


37.0 


10.12 ± 0.04 


1991.73 ± 0.01 


1.36 ± 0.01 


0.00 ± 0.02 


0.76 ± 0.14 


1.91 ± 0.04 




22.0 


9.72 ±0.07 


1991.99 ±0.02 


1.55 ±0.03 


0.26 ±0.03 








14.5 


7.42 ±0.02 


1992.65 ±0.01 


2.12 ±0.01 


0.91 ±0.01 








8.0 


4.52 ±0.03 


1992.88 ±0.01 


1.51 ±0.01 


1.15 ±0.01 








4.8 


2.15 ± 0.01 


1993.31 ± 0.01 


1.82 ± 0.01 


1.58 ± 0.01 






F 


37.0 


5.45 ±0.04 


1992.97 ±0.01 


0.99 ±0.01 


0.00 ± 0.01 








22.0 


4.08 ± 0.06 


1993.10 ± 0.02 


0.85 ± 0.02 


0.13 ± 0.02 






G 


37.0 


1.83 ±0.03 


1993.93 ±0.01 


0.80 ±0.01 


0.00 ± 0.02 








22.0 


2.46 ± 0.03 


1993.93 ± 0.01 


0.86 ± 0.01 


-0.01 ± 0.02 






H 


37.0 


4.36 ± 0.03 


1995.57 ± 0.01 


2.17 ± 0.03 


0.00 ± 0.02 


0.13 ± 0.07 


1.46 ± 0.03 




22.0 


5.30 ± 0.03 


1995.52 ± 0.01 


1.99 ± 0.01 


0.05 ± 0.02 








14.5 


4.08 ±0.02 


1995.69 ±0.01 


2.43 ±0.01 


0.12 ±0.02 








8.0 


4.17 ± 0.02 


1995.61 ± 0.03 


4.06 ±0.05 


0.04 ± 0.04 








4.8 


3.51 ±0.02 


1995.97 ±0.02 


3.96 ±0.06 


0.40 ± 0.03 






I 


22.0 


1.27 ±0.02 


1997.18 ±0.01 


1.13 ±0.01 


0.00 ± 0.03 


-0.32 ±0.19 






14.5 


1.81 ±0.08 


1998.06 ± 0.07 


1.40 ±0.04 


0.88 ± 0.08 








8.0 


1.46 ±0.03 


1998.69 ±0.01 


1.53 ±0.01 


1.51 ±0.03 








4.8 


2.41 ±0.01 


1999.01 ±0.02 


1.92 ±0.01 


1.83 ±0.02 






J 


37.0 


3.23 ±0.04 


1998.62 ±0.01 


1.96 ±0.01 


0.00 ± 0.02 


0.17 ±0.27 


0.40 ±0.01 




22.0 


3.65 ±0.03 


1998.82 ±0.01 


1.70 ±0.02 


0.20 ±0.02 








14.5 


2.70 ±0.02 


1999.32 ±0.01 


1.26 ±0.01 


0.70 ± 0.02 






K 


37.0 


5.51 ±0.05 


2000.97 ± 0.02 


3.35 ±0.05 


0.00 ± 0.03 


-0.17 ±0.01 






22.0 


6.26 ±0.06 


2001.27 ±0.03 


3.01 ±0.09 


0.30 ± 0.05 








14.5 


6.46 ±0.01 


2001.24 ±0.02 


2.51 ±0.01 


0.27 ±0.02 








8.0 


7.25 ±0.02 


2001.16 ±0.01 


3.00 ±0.01 


0.19 ±0.02 








4.8 


7.93 ±0.01 


2001.38 ±0.01 


2.49 ±0.01 


0.41 ±0.02 







Table 2. Comparison of frequency-dependent core shifts measured from VLBI observations and calculated from the frequency-dependent 
time lags for 3C 345. 



Core shift from VLBI 


Core shift from time lags 


Epoch 


Jet 


(mas) 


(mas) 




comp. 


0.05 ±0.13 (10.7 - 5 GHz) 


0.06 ± 0.02 (4.8 - 8 GHz) 


1982 


C2 


0.111 ± 0.007 (22.2 - 10.7 GHz) 


0.12 ±0.02 (22.2 - 8 GHz) 


1983.5 


C4 


0.05 ±0.03 (5.0 - 8.4 GHz) 


0.09 ±0.01 (4.8- 8 GHz) 


1993.8 


C7 


0.21 ±0.06 (8.4 - 22.2 GHz) 


0.19 ±0.02 (8.0 - 22.2 GHz) 


1992.8 


C7 


0.33 ±0.10 (5.0 - 22.2 GHz) 


0.28 ± 0.04 (4.8 - 22.2 GHz) 


1993.7 


C7 
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22.2 GHz), and Ar = 0.33 ± 0.10 mas (5-22.2 GHz) in 
~ 1993 (see Table [2j. We calculated core shifts using our 
measured time lags for the flare in 1992: 0.43 ± 0.01 yrs 
(4.8-8.0 GHz), 0.89 ± 0.03 yrs (8.0-22.2 GHz), and 1.32 ± 
0.03 yrs (4.8-22.0 GHz). The speed of the jet component 
C7 ejected during the 19 92 flare has been estimat ed to 
be 0.208 ± 0.025 m as/vr dCaproni fc Abraham! \2004 ) and 
0.30 ± 0.16 mas/yr l|Ros et al J 120001 ). Taking this speed for 
the newly ejected jet component C7 and using the mea- 
sured time lags, we can calculate the frequency-dependent 
core shifts, shown in Table [2] 

Table[2]clearly shows that the core shifts measured from 
the VLBI observations and calculated using our frequency- 
dependent time lags and the measured jet-component speeds 
coincide very well. This provides a direct test of the proposed 
method, and suggests that it can be used to reliably calculate 
core shifts from total flux-density light curves. Moreover, 
this provides direct evidence that the jet component speed 
is the actual physical speed of a knot or a shock, rather than 
a pattern speed. 
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Figure 4. 3C 345: Plot of time lag measurements versus fre- 
quency, using 37 GHz as the reference frequency. The solid curve 
shows the best-fit line At = a(i^ _1 /*'-) _ 37.o(~ 1 / fc '')) with a = 
5.47±0.04 and k r = 1.91±0.04. 



5 RESULTS FOR THE QUASAR 3C 345 

Looking at Eq. (3), we expect a plot of the time lags versus 
frequency to have the form 

At(u) bs = Constant * (i>a^^ kT — ^fT r )- (12) 

The plots of time lags versus frequency for 3C 345 show that 
these do, indeed, follow such a power law, enabling us to fit 
the time lags and derive k r values from the fits. Figure [3] 
shows the fit for flare E as an example. The estimated k r 
value from this plot is 1.91L0.04. We used the highest fre- 
quency, 37 GHz, as the reference frequency for all measured 
time lags. Measuring k r can tell us about the jet geometry, 
since the distance from the VLBI core to the base of the 
jet scales with k r as r core ~ u^ 1 ^ 1 ' (Eq. [5]). The k r values 
contain information about the distributions of the magnetic 
field and electron number density, since k r depends on the 
indices m and n (Eq. [SJ, which indicate how the electron 
number density and magnetic field decrease along the jet 
(Eq. [8}. The calculated k r values for each flare of 3C 345 
are listed in the last column of Table Q] The range of our 
k r values enc ompasses t he av erage value k r = 0.96 derived 
for 3C 345 bv lLobanovl (I1998T ). We can also use k r to calcu- 
late the core-region magnetic fields using Eqs.m fTITl and llll 
Our k r values imply magnetic fields in the core of 3C 345 
~ 0.1 G. 

The measured k r values for 3C 345 suggest a significant 
time evolution. The time delays versus frequency are evolv- 
ing with time. Figure [5] shows time lags for individual flares. 
The maximum time lag is changing between 0.4 yrs and 1.8 
yrs. Figure |S] shows that the k r values are almost at the 
same level in the period from 1982 to 1992, about 1.8, then 
begin to decrease to k r ~ 0.4 in 1999. In this same period, 
1982-1990, 3C 345 experienced a major flare, reaching flux- 
density levels of about 19 Jy (Fig. [2j . The lower k r values 
seem to correspond to flares with less dramatic amplitudes, 
suggesting a possible connection between k r and the flux- 
density level. Figure shows k r versus the flux at 14.5 GHz; 
there is a clear tendency for brighter flares to have higher k r 
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Figure 5. 3C 345: Time delay of individual events as functions 
of frequency. 

values, reaching a kind of saturation level of about k r ~ 1.8 
for fluxes higher than 6 Jy at 14.5 GHz. 



6 DISCUSSION AND CONCLUSION 

We have introduced a new method for calculating core shifts 
from time lags between the maxima of single-dish flares at 
different frequencies. The method assumes that the observed 
velocity of jet components is the actual speed of matter in 
the jet. 

Applying the method to the 3C 345 VLBI jet and inte- 
grated light curves from 4.8 to 37 GHz, we found that the 
directly measured frequency-dependent core shifts and the 
core shifts calculated with our new method agree very well, 
supporting the validity of the method and our assumptions. 
In particular, this provides direct evidence that the observed 
component speeds in the VLBI jet represent the actual phys- 
ical speeds of these components, rather than the pattern 
speed of a shock. This technique should also be checked us- 
ing more sources, which we plan to do in a future study. 



Frequency- dependent core shifts in AGN jets. 

Table 3. Calculated cquipartition magnetic fields and distances between the radio core and the base of the jet. 
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v Ar fl rv r core (v) B lpc B core (u) 

[GHz] [mas] [pc-GHz] [pc] [G] [G] 

1641+399 (3C 345): Flare H, k r = 1.46 ±0.03, /3 app = 0.37 ± 0.03 1 

4.8/37 0.15 ±0.02 3.79 ± 0.44 6.79 ±0.78 0.46 ± 0.08 0.07 ± 0.01 

8/37 0.01 ±0.01 0.59 ±0.55 1.06 ± 0.99 0.11 ±0.16 0.11 ±0.18 

14.5/37 0.043 ±0.007 3.76 ± 0.61 6.73 ±1.10 0.46 ±0.11 0.07 ± 0.02 

22/37 -0.019 ±0.007 3.50 ± 1.29 6.26 ± 2.31 0.43 ± 0.24 0.07 ± 0.05 



Col.(l): frequency; Col. (2): calculated core shift from frequency-dependent time 



delays; Col. (3): Q. r u\ Col. (4): distance from radio core to the base of the jet; Col. (5): 
cquipartition magne tic field at 1 pc; C ol. (6): equipartition magnetic field in the core. 
References: 1 iKellermann et al.l ll2004h . 
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Figure 6. 3C 345: Time evolution of k r values. 
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Figure 7. 3C 345: correlation between k r values and amplitude 
of the flares at 14.5 GHz. The picture shows that the k r values 
become higher with higher flux levels. 



We have used our method to derive k r values for 3C 345. 
We find clear evidence for variability of k r , with the mea- 
sured values ranging from 0.4 to 1.9. We find evidence that 
k r increases with the core flux level, reaching saturation at 
a value of ~ 1.8 above core fluxes of about 6 Jy. In principle, 
time evolution of k r could come about due to changes in the 



core spectral index, magnetic-field distribution, or electron 
number-density distribution, since k r depends on a, m, and 
n (see formula [5]). There is no obvious relationship between 
the value of k r and the core spectral index (Table [TJ. Un- 
fortunately, we cannot directly separate the m and n values 
from the k r equation (the only exception is if k r is close 
to unity, i.e., close to equipartition, in which case it is rea- 
sonable to infer m = 1 and n = 2). Therefore, we cannot 
unambiguously prove that k r variations are due, for exam- 
ple, purely to variations in a, m or n; it is likely that all three 
parameters contribute to time variability of k r . Time vari- 
ability of k r values implies that frequency-dependent core 
shifts are changing with time, suggesting that in order to 
match the astrometric catalogues it is necessary to have si- 
multaneous multi-frequency observations. 

Using our k r values we can estimate the distance of 
the radio core from the base of the jet, the equiparti- 
tion magnetic field at 1 pc distance, and the equi partition 
magne tic field in the core (|Lobanovl (1 19981 ) and iHirotanH 
(2005)). Since not all outbursts are in the equipartition 
regime as was shown in the previous sections, we have se- 
lected for magnetic fields calculations one outburst H with 
k r value from a Table [1] closest to unity (and therefore 
to equipartition regime). We have used intrin sic jet half- 
opening ang le 9 = 0.5° Jjorstad et all [iooBh . jet view- 



20051) . Doppler factor 



ing ang le ip = 2.7° dJorstad et al.l _ 
5 — 7.8 (|Hovatta et al.ll2009l ). and luminosity distance Dl = 
3473 Mpc. Table |3] shows calculated frequency-dependent 
core shifts, distances between the radio core and the base 
of the jet, magnetic field at 1 pc, and magnetic field in the 
core for individual pairs of frequencies. We have not used 
8 GHz data in the analysis, since the frequency-dependent 
time lag for the 8/37 GHz pair of frequencies is not very 
reliable. In error analysis we took into account only the 
errors in time lags measurements and apparent speeds of 
jet components. The averaged magnetic field in the core of 
3C 345 is B CO re = 0.07 ± 0.02 G and magnetic field at 1 pc 
is B lpc = 0.45 ± 0.09 G. 

By calculating the frequency-dependent core shifts from 
the frequency-dependent time delays measured for inte- 
grated light curves, we can study the long-term evolution 
of the core shifts and calculate the core shifts for any par- 
ticular time when long-term radio light-curves are available. 
Since the total flux-density light curves covering more than 
30 years are available for dozens of radio sources, this makes 
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it possible, in principle, to calculate the core-shift evolution 
over more than 30 years without constructing and aligning 
VLBI maps at multiple frequencies. The only input needed 
from direct VLBI observations is the apparent speeds of jet 
components (which can be measured from a series of obser- 
vations at a single frequency). 
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